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Detection of differential rotation in V^ Cap with profile analysis^ 
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Abstract. We report detection of differential rotation on the F5 dwarf t/j Cap using line profile analysis. The 
Fourier transform of both Pel A5775 and Sil A5772 are used to obtain a projected rotational velocity of v sini — 
42 ± 1 km s~ ^ . Modelling of the Fourier transformed profiles shows that the combined effects of equatorial velocity, 
inclination and differential rotation dominate the line profile while limb darkening and turbulence velocities have 
only minor effects. Rigid rotation is shown to be inconsistent with the measured profiles. Modelling the line profiles 
analogous to solar differential rotation we find a differential rotation parameter of a = 0.15 ± 0.1 (15 ± 10%) 
comparable to the solar case. To our knowledge this is the first successful measurement of differential rotation 
through line profile analysis. 
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1. Introduction 



Accor ding to Llie staudaiJ paiadigui of sLuUai activity 
differential rotation is a central ingredient for the mag- 
netic dynamo presumed to underlie all activity phenom- 
ena. In a turbulent convection zone differential rotation is 
expected to be a central element of stellar activity. The 
interaction of rotation and convection naturally produces 
deviations from rigid rotation. One example for this ef- 
fect is the Sun, who se equator rotates ~ 20% faster than 
higher latitudes (e.g. Lang 1992 ). Direct predictions of the 
dependence of differential rotation on ro tational velocity 
and spectral type have been made (e.g. Belvedere et al 



198C ; Kitchatinov fc Riidiger 1999 ). Observations suggest 



a strong correlation between rotation and activity in late- 
type stars with many pieces of evidence summarized as 
"rotation-activity c onnection" ( Hempclmann et al. 1995 ; 
Messina et al. 200X1 ). 

However, measurements of differential rotation among 
stars are either complicated or time-consuming or both. 
One method to determine differential rotation is to exam- 
ine spectral line shapes. Subtle differences exist between 
the line profile of a rigidly rotating star and that of a dif- 
ferentially rotating star, because the underlying flow fields 
are different. Although significant differences between the 
two cases are predicted, all previous attempts to measure 
differential rotation through line profiles alone remained 
unsuccessful and yielded results consistent with rigid ro- 
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* Based on observations collected at the European Southern 
Observatory, La Silla (65.L-0101) 



tation (e.g. [Gray 1982 ; Dravins et al. 1990; Groot et 
al. 199q). At any rate, for a successful measurement one 



needs high spectral resolution and high signal to noise. 
With the Very Long Camera (VLC) of the Coude Echelle 
Spectrograph (CES) at ESO's 3.6m telescope an instru- 
mental setup satisfying these requirements is available. In 
this letter we report the detection of differential rotation 
in the F5 dwarf ip Cap, to our knowledge the first suc- 
cessful measurement of differential rotation through line 
profiles. 



2. Data 



4.13, 



V-Cap (H P 197692, F5V, K nag 

40kms^^; Ucsugi et al. (1982 )) is the fastest rotator in 



a set of solar-like stars in the solar neighbourhood we 
observed on October 13, 2000 at ESO's 3.6m telescope 
(La Silla). The spectral resolution achieved with the CES 
/ VLC setup was R ^ 235 000 (- 1.28kms"i). Three 
consecutive exposures of ^Cap of 270 s each covering 
the wavelength range between 5770 - 5810 A were taken. 
The signal to noise ratio of the extracted spectrum is 
S/N ~ 800 per pixel with proper fiat fielding and removal 
of interference pattern (c.f. Kiirster 2001, Sect. 6). 

We concentrate our analysis on Fourier transforms of 
wefi isolated lines (FeiA5775 and SiiA5772), for which 
continuum placement because of line blending is not a 
problem even for large v sin? values. The signal to noise 
ratio of our data is high enough to show the crucial fea- 
tures in the Fourier transform already from a single line. 
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that Doppler broadenings can be treated as convolutions; 
and that the "intrinsic" profiles are identical over the ap- 
parent stellar disk. For fast rotators {v sini > 30 km s^^) 
and stationary atmospheres rotational broadening domi- 
nates and no complications occur with this assumption. 

Interpreting the observed data profile D{X) as a a mul- 
tiple convolution (denoted by *) between the intrinsic pro- 
file F{X) (including microturbulence broadening), the ro- 
tational broadening profile G{X), the instrumental profile 
/(A) and the macroturbulence profile M{X), D{X) can be 
written as 




\. 



D{X) = F{X) * G(A) * /(A) * M{X). 



(1) 



We calculated F{X) using the packages Atlas9 
( iKurucz 1979| |1993| ) and Bht ( [Baschck el al. 1966|) 
Atomic damping coefficients obtained from VALD (Kupka 



et al. 1999; Kurucz 1994) were included in the profile cal- 



Fig. 1. (a) FeiA5775 (right) and SiiA5772 (left) lines of 
tp Cap (upper line) and t Psc (lower line; shifted along y- 
axis) . (b) Original and modified line profiles of ^ Cap (see 
text). 

Still, in FeiA5775 a small blend occurs from a neigh- 
bouring set of weak lines ^ 0.8 A apart, the influence of 
which on differential rotation determination is not clear. 
During our observing run a spectrum of the slow rota- 
tor lPsc (HD 222 368; F7V) was taken. In Fig.^ both 
spectra are shown, with the profile of t Psc shifted accord- 
ing to its relative radial velocity. Deblending of Fei A5775 
was accomplished by, first, removing the line in the spec- 
trum of iPsc "by hand", broadening this template with 
Wrot = 40kms"^ and subtracting it from the V'Cap spec- 
trum. In Fig. Hp the recorded spectrum as well as the 
modified spectrum of ip Cap are shown. Obviously, the 
SiiA5772 equivalent widths of ipCap and iPsc are the 
same, since it has been successfully removed (cf. FigMj). 
The Fe l A5775 line appears symmetric, thus we think that 
deblending was successful. 



3. Method of Analysis 

Absorption line profiles are influenced by a number of dif- 
ferent effects. An absorption line profile at any point on 
the star is determined by temperature, gravity, element 
abundances and atomic parameters. This "intrinsic" pro- 
file is Doppler broadened by velocity fields. Many efforts 
have been undertaken to distinguish these velocity fields, 
one of them being the stellar rotation (see Gray 1988, and 



references therein). In addition to the projected rotational 
velocity of the star, radial-tangential macro- and isotropic 
microturbulence (denoted with Crt and ^ resp.) turned 
out to be a reasonable parameterization of stellar atmo- 
spheric velocity fields. These parameterizations assume 



culations; we chose solar metallicity and for the velocity 
dispersion of the opacity distribution function in Atlas9 
we used 1.0km s~^. The temperature has been set to 
Tcff = 6500 K dBlackwell fc Lynas-Gray 1998|). To cal- 



culate G(A) a modified version of a package developed 
and described by Townsend (1997) is used. The surface 
integration is carried out over about 25 500 visible surface 
elements. The adopted limb-darkening law is given by 



I{cos9) = Io{l — e + ecosl 



(2) 



with denoting the angle between the surface element 
normal and the observer's line of sight and e the limb 
darkening coeSicient. We parametrize the differential ro- 
tation law through 



uj{l) = ujQ — uji sin I 



(3) 



with / being the latitude. Specifically, differential rota- 
tion is expressed in terms of q = loi/uq. The differen- 
tial rotation law (Eq.|^) is adopted from the solar case, 
where aQ = 0.20; a more general approach would ex- 
pand uj{l) in terms of orthogonal polynomials. The instru- 
mental profile /(A) has been determined from the shapes 
of calibration lamp lines, a correction for thermal broad- 
ening according to the lamp temperature of 120°C was 
applied. The m acroturbulen ce broadening function M{X) 
is adopted from [Gray (1992i p. 407). 

The contributions of different velocity fields are dif- 
ficult to separate especially in the wavelength domain. 
In the Fourier domain convolutions become multiplica- 
tions which are much easier to handle; the advantages of 
Fourier domain are discussed in detail by Gray (1992 ). In 
the Fourier domain Eq. H^ becomes 



d{cr) = /(f^) • 5(0-) • «(cr) • m{a). 



(4) 



The Fourier frequency a is expressed in cy- 
cles / (km s^^) — s km^^. Noise in Fourier domain 
can be expressed as Sa = S^AXN"! ( Gray 1992 ), 
with Sx being data noise and AA step size in the 
wavelength domain, N is the number of data points. 
The Fourier transform of a real, symmetric profile yields 
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a real function in the Fourier domain. Investigations 
of imaginary transforms of asymmetric lines show that 
these are no reliable tracers of profile properties (Gray 
1980|)T]We created a symmetric profile by mirroring the 
absorption profile at it's center; the achieved Fourier 
transform turned out to be stable against small variations 
of the center position. Rotational broadening, which 
dominates the line broadening in the case of ip Cap, yields 
symmetric profiles. Asymmetries due to convection are 
believed to be of the order v < lkms~^. In Fourier 
space this is cr > Iskm"'^ while our analysis focuses on 
a < 0.04skm^^. In the case of FeiA5775 an asymmetry 
may occur due to imperfect deblending, for SiiA5772 no 
mechanism should contribute asymmetries of this order. 

The Fourier transformed observed line profile was com- 
pared with a Fourier transformed model profile via a x^ 
test. Within our adopted modelling approach the following 
six fit parameters determine the model profiles: rotational 
velocity (u), inclination angle (i), differential rotation (a), 
limb darkening (e), macro- (Crt) and microturbulence (^). 
The x^-calculations have always been carried out directly 
on d{<7). No further corrections have been applied to the 



data to avoid complications with amplified noise. Bruning 
(1984|)]pointed out that the convolution method induces 
systematic errors due to incorrectly estimated line depths 
especially in slow rotators. Our study is not focused on re- 
producing equivalent widths, which can be tuned e.g. by 
element abundances. The main goal is to reproduce the 
line shapes. Although we do not examine slow rotators 
in this study, normalized transformed profiles were used, 
that is, Fourier transformed profiles have been scaled to 
unity at a = 0. 

4. Results and Discussion 

After a rough determination of v sini a grid calculation 
in the six fit parameters was carried out. As to our model 
parameterization we note that this description contains 
parameters whose physics is poorly understood, i.e., the 
micro- and macroturbulence. We thus use x^ only as a 
relative parameter and not as an absolute criterion as to 
whether a model is acceptable or not; Xn denotes the good- 
ness of fit relative to the best fit, for which x^ = 1 is set. 
For the actual grid we chose the parameter ranges of e (0.4 
- 0.8), Crt (4.0 - 7.0) kms^^ and f (1.0 - 2.5) kms^^. 
The rotation velocity was limited to w < 120 km s^^ 
(i > 20°). 

Goodness of fit calculations have been carried out on 
FeiA5775 and Sii A5772; both Hues show identical results. 
In the following we show only the results for FeiA5775. 

In Fig. II we plot the resulting Xn as a function of the 
differential rotation parameter a; for each fixed value of a, 
we plot the Xn values obtained by varying all other model 
parameters. Clearly, a well defined lower envelope curve 
exists with a minimum at a = 0.20. The best fits for dif- 
ferential rotation (here a = 0.20; Xn — 1j Model 1) and 
rigid rotation {a — 0.0, Xn = '^•9, Model 2) are marked 
by full squares; all parameters are shown in Tab.^. 
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Fig. 2. Normalized x^ vs. a for calculated models with 
various values of «,«,€, Crt and £, (see text); each cross 
represents one set of parameters. Best fits of all models 
(1) and for rigid rotation (2) are denoted by full squares. 

Table 1. Parameters of best fitting models for rigid and 
differential rotation 



Model 



Crt C 



^ 60° 0.00 0.8 7.0^^ 2.5^ 7.9 




Fig. 3. The Fourier transformed data indicated by error 
bars. Best fits for differential rotation {a — 0.20; full 
line. Model 1) and rigid rotation {a — 0.0; dashed line, 
Model 2) are shown. See Tab. [I] for parameters. 

In Fig.y we plot the corresponding profiles in the 
Fourier domain; obviously. Model 1 provides a much bet- 
ter fit than Model 2, which is not an adequate description 
of our data. While no statistical uncertainties can be given 
on our Xn-values, Fig. g shows xfr'^^lnss of 7.9 to be clearly 
unacceptable. 

Inspection of the calculated Xn'S^'^^ shows that the pa- 
rameters e, Crt and ^ produce only second order effects; 
the fits are driven by the chosen values of v, i and a, 
and we cannot determine all model parameter s indepen- 
dently Therefore, we fixed the values e = 0.6 ( [Carbon fc 
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and PZTel (Barnes et al. 200C). As direct predictions 
for a F5 dwarf have not been calculated by Kitchatinov 



ji.^. 
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Fig. 4. Xn of fits to Fe I A5775 with 41 .3 km s-1 


< V sini < 


43.3kms-\ e = 0.6, Crt = 6.0kms-i 


and ^ — 


l.Okms^^ in the a-i (v) plane: Xn < 5(»); 


5 < Xn < 


10(o);x2>10(.). 





Gingerich 1969), Crt = 6.0kms~^ and ^ = l.Okms"^ 
( pray 198^ , and references therein). We emphasize that 
our subsequent findings do not require these specific pa- 
rameter settings, and that the influences of e, Crt and £, are 
too small to revoke the effects oi a, v and i. For the given 
parameters the best fit value for v sini is 42.3 km s~^. We 
estimate the total systematic errors to be ~ lkms~^. 
Models outside the range 42.3 ± 1 kms"""^ all have xi > 5. 
Theoretical investigations of the line profile behaviour 
in Fourier space show that with differential rotation a 7^ 
the inclination i becomes important (e.g. Bruning 1981). 



For comparable values of v sin i the best fit values of a and 
i are correlated. In Fig. |j we consider calculated models 
with fixed e, Crt, C, ''J sini between 41.3 and 43.3kms^^ 
and varying a, v and i. Three groups of models are distin- 
guished; xl < 5(»), 5 <xl < 10(0) and xl > 10(-)- 
A well defined area of reliable fits in the a-i plot emerges. 
All fits with Xn < 10 show a > 0.0. Although we used no 
absolute criterion on wheter a model is acceptable, the es- 
timation of systematic errors and a comparison with Fig. ^ 
shows that the threshold Xn < 10 is a rather high choice. 

A variety of combinations of equatorial velocity and in- 
clination seems possible for a > 0. A smaller differential 
effect is preferred for faster rotating models. For extremely 
high values of v > 200kms~^ rigid rotation might be ar- 
gued, but such velocities seem unlikely given the age of 
tjj Cap of '^ 2 10^ a (Lachaumc et al. 1999). Even the min- 
imum rotational velocity of 42kms~^ seems uncommonly 
high. From Fig. we derive a = 0.15 ± 0.1. The identical 
result was found for SiiA5772, where the same procedure 
was applied. 

To summarize differential rotation has been estab- 
lished for the rapid rotator -0 Cap independently from two 
absorption fine profiles. While '0Cap rotates at least 20 
times faster than the Sun, it's differential rotation is com- 
parable to the solar value, but not with the differential 
rotation patterns deter mined from Doppler images of the 
fast rotators ABDor ( Donati fc Collier Cameron 1997 ) 



& Riidiger (199£) and the rotation period of ■0 Cap is 
only poorly determined, the consistency of our result 
and the model is not clear. Assuming B = 1.2 Rq and 
i — 90°, we find for -0 Cap the rotation law uj{l) = 
4.38 — 0.66 sin^ I rad/d; V" Cap does not rotate like a rigid 
body as suggested for AB Dor and PZ Tel. To compare 
theory and observations more detailed predictions espe- 
cially on spectral type dependence and a greater sample 
of direct observations are needed. In particular verification 
of the differential rotation results of line profile analysis 
and Doppler imaging for the same star will be instructive. 
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